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Abstract 

An analytic expression is derived for the distribution G(R) of the 
end-to-end distance R of semiflexible polymers in external potentials 
to elucidate the effect of confinement on the mechanical and statisti- 
cal properties of biomolecules. For parabolic confinement the result is 
exact whereas for realistic potentials a self-consistent ansatz is devel- 
oped, so that G(R) is given explicitly even for hard wall confinement. 
The theoretical result is in excellent quantitative agreement with fluo- 
rescence microscopy data for actin filaments confined in rectangularly 
shaped microchannels. This allows an unambiguous determination 
of persistence length Lp and the dependence of statistical properties 
such as Odijk's deflection length A on the channel width D. It is 
shown that neglecting the effect of confinement leads to a significant 
overestimation of bending rigidities for filaments. 

Bending of semiflexible macromolecules such as actin filaments play a 
crucial role for the mechanical properties of cells pp. In general, the motion 
of biopolymers such as DNA or actin filaments takes place in gels or cyto- 
plasm where sterical constraints forces a single molecule to bend in addition 
to the bending constantly induced by thermal motion. Advances in microflu- 
idic techniques and in the direct visualization of actin filaments [21 El HI El E] 
make it nowadays possible to study experimentally the interplay of ther- 
mal fluctuations and confinement in controlled environments. Although the 
dynamics of single biomolecules can be measured for almost 20 years by flu- 
orescence microscopy in simple confining geometries [TJ [HI El CD] as well as 
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in microfabricated porous arrays [11, 12J, the theoretical understanding of 
such constrained thermal motion is still hindered by the fundamental statis- 
tical problem to treat thermal activated undulations of macromolecules and 
steric repulsion from obstacles on the same analytical level. Whereas ther- 
mal bending modes can be treated well in Fourier space and steric constraints 
in real space, respectively, the combined interaction of thermal bending in 
confinement is not solved yet in a satisfactory way. 

In order to characterize thermal fluctuations of polymers of length L one 
may calculate the tangent-tangent correlation function or the radial distribu- 
tion function, i.e., the probability distribution G(f) for the distance vector r 
of the two filament ends [6]. Because end-to-end distances can be measured 
easily by labelling the ends of the macromolecule the latter served in the 
past as an important experimental tool to elucidate the physical properties 
of biopolymers including DNA in nanochannels [T3J. However, estimated 
values for bending rigidity k and persistence length Lp = K,/(kp,T), for in- 
stance, depend on the details of confinement. An expression for G(r) for 
freely fluctuating single, semiflexible polymers can be found in Ref. [Hj, 
which is unfortunately not applicable for confined and interacting filaments. 
The effect of confinement on the radial distribution function G(r) is also rel- 
evant for the microrheology of actin solutions where entanglement plays an 
important role |15j . Thus, an analytical result for the dependence of G(r) on 
confinement would provide an important tool to analyze experimental data in 
microfluidic devices in order to determine unambiguously the physical prop- 
erties of biomolecules. The statistical mechanics of unconstraint fluctuations 
of semiflexible chain molecules (L < Lp) is well understood. A remarkable 
successful description is the worm-like chain model based on a fluctuating, 
elastic string |16| . However, not even the case of parabolic confinement was 
solved yet beyond tangent correlations functions [T7j. Here, we present an 
analytic solution for the end-to-end distribution function for strong confine- 
ment with D « L which is in quantitative agreement with experimental 
findings measured in Refs. p3, El EE], for instance. Moreover, a self-consistent 
ansatz allows the unambiguous mapping of non-parabolic confinements like 
hard walls or even van-der-Waals potentials on a quadratic Hamiltonian [TH] , 
so that our analytic result is even applicable on realistic substrate-filament 
interactions. This self-consistent ansatz was already successful in predicting 
the fluctuation spectrum of bounded membranes in excellent agreement with 
experimental data [18J as well as for thickness dependence of wetting layers 
|19j . Here, we present the first application on one-dimensional fluctuating 
filaments and compare our analytical result with experimental data for actin 
filaments [6]. 

Actin filaments have diameters of ~ 8nm, lengths of L = 11 — 13/im and a 
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Figure 1: Microfluidic devices fabricated by soft photo-lithography makes 
it possible to observe thermal fluctuations of actin filaments in confining 
microchannels by fluorescence microscopy (image scanned from Ref. 
A theoretical description can be based on a Monge-parametrization of the 
perpendicular deviation f(z) of the filament from the z-axis along the channel 
instead of its arc length s. Despite the limited spatial resolution of the 
image the end-to-end distance R can be measured accurately which makes 
its distribution G(R) an ideal quantity to determine physical parameters of 
the filament. G(R) is given analytically in Eq. ([3]) in very good agreement 
with experimental data shown in Fig. [HJ 

persistence length of about Lp « 15/mi yielding thermal fluctuations which 
can be observed by optical fluorescence microscopy [UJ. Thus, in addition to 
the importance of these thermal fluctuations for the biological functioning 
of cells, actin filaments can be used as an experimentally accessible model 
system for semiflexible polymers where L < Lp. The statistical properties 
of confined filaments depend crucially on the ratio L/X of the total length 
L of the filament to Odijks deflection length A which equals approximately 
the number of bendings caused by the walls. Because the channel width D 
in the experiments described in Ref. [UJ is much smaller than the length of 
the macromolecules we focus here on the strong confinement limit D « 
L. An actin filament of length L is modelled as a differentiable curve in 
space parametrized by the position vector r(s) and the arc length s G [0, L] 
(see Fig. [T]). Because polymers in microchannels are elongated as shown 
in the fluorescence microscopy image, one may ignore loops or overhangs of 
f(s). Applying a Monge parametrisation one can write the thermally induced 
deviation f(z) = (f x (z), f y {z)) from the z-axis e z parallel to the channel walls 




Then, the Hamiltonian of 



an elastic string with bending rigidity k can be written as 




(1) 
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with the normalized tangent vector t(s) = %. = (d z f, lj / \J 1 + (d z f) 2 , the 
abbreviation d z = J^, and the projection L z of the contour length L on the 
z-axis of the channel. Notice, that L z (f) itself is a function of / given by 

the implicit relation L = f^*^ dz y 1 + (d z f(z)) 2 at fixed contour length 
L. The channel walls are assumed to be purely repulsive with a quadratic 
cross-section and width D. To model the confinement we add an energy term 
7~Lpot = \ Jo z dz U(f(z)) where the potential U describes the interaction of 
the filament with the substrate material. 

In order to avoid non-Gaussian path integrals in the beginning, we start 
with a parabolic potential U(f) = -f/ 2 perpendicular to the channel axis 
e z . In a second step we map self-consistently non-Gaussian potentials on 
an effective parabolic potential strength E and derive an explicit expres- 
sion for the dependence of E on the channel width D. Accordingly, we 
expand also TLbend up to second order in / because the fluctuations f x ,y{ z ) 
are expected to be small at room and even physiological temperatures, yield- 
ing H = 7Ybcnd + Upot = \ Jq Z dz K,(d z f) 2 + Ef 2 . Thus, the statistical 
properties of the actin filament are determined by only two characteristic 
length scales, the persistence length Lp = «/(/csT) and the deflection length 

A = £|| = introduced by Odijk [32]. Whereas L P is a quantitative mea- 
sure for the thermal flexibility of the filament, A is the characteristic length 
of bendings forced by the confinement. The Gaussian approximation is ad- 
ditionally justified by the channel confinement which suppresses effectively 
large fluctuations on length scale larger than A. Notice, that the inextensi- 
bility constraint \i\ — 1 is violated in the Gaussian approximation, although 
automatically fulfilled by using the arc length s in t = |j. 

We exploit open end boundary conditions and mirror the polymer at the 
origin in order to be able to apply periodic boundary conditions for this aug- 
mented configuration. The mirroring of the polymer suppresses sine-modes 
but doubles the cosine-modes by allowing wavevectors jk (k G N) instead of 
^■k as for usual periodic boundary conditions. By doing so, the end and the 
midpoint of the periodic string, i.e. the two endpoints of the original filament 
can fluctuate freely and are not constrained on the same excursion, which 
is crucial for the following analytic calculation. This boundary condition is 
equivalent to the one used by Wilhelm and Frey for unconfined polymers [2] , 
so that we recover their result for E = as one can see below. Then, f(z) 
can be expressed in a discrete Fourier sum f x ,y{z) = Ylh=i /W cos {lT^ z 
with integer k. The end-to-end vector R = ds t(s) 
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reads in Gaussian approximation in terms of the Fourier amplitudes 



*=«.-E( (1 ^7™) +ocftt)1) ( 2 ) 



fc=i 



yielding the end-to-end distribution function G(r) = (^5 — Rj^ = 6(r x )8(r y )Go(r) 

where the angle brackets indicate an ensemble average determined by the 
Hamiltonian JT]) in Gaussian approximation. Applying Fourier expansion of 

5 ( f — R ) as well as residuum theorem one finds 



0[ ) LM^ w A A F(k')-F(k) 

k=l fe'6N\{fc} V ' K ' 



^Kf-i.^WlAM) (3 ) 
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with the dimensionless length ratios 

K = ~ and F(fc) = ~f-n 2 k 2 + (4) 

An essential step in the derivation of this analytic result is the splitting of 
the product 

n F{k ' ] = n k ' 2 n 1 + ^ (5) 

k'^k y ' y ' k'^k k'^k 1 k' 2 k 2 

into two factors. The first product equals (— l) fc_1 whereas the second one 
can be rewritten in terms of the functions p n {x) — Yl ( 1 + Ir ) with P2(x) = 

.7=1 V ^ 1 



sin ttW—x 



-. Because Pi{k^) is a numerical factor independent on k, it can be 
put into the normalization constant M which is determined by J d 3 r G(r) = 
1. 

In Fig. El we illustrate the dependence of Gq(t) on the potential strength 
E as well as on the bending rigidity k, i.e., on A and Lp, respectively. We 
observe that with increasing E the most probable end-to-end distance r/L, 
i.e., the position of the peak in Gq(t) is shifted to higher values and that the 
peak gets sharper. This is obvious since strong confinement supports rather 
straight conformations and short end-to-end distances get more unlikely. In 
the limit of vanishing confinement E — > one finds a diverging deflection 
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Figure 2: Radial distribution function Gq(R) for filaments of persistence 
length Lp as function of the confining potential strength E and the contour 
length L. 



length A — > oo and a vanishing critical wavevector k c — > 0. In this limit, Eq. 
([3]) is identical to the result of Wilhelm and Frey derived for freely fluctuating 
semiflexible polymers [H] where large wavevectors k do not contribute signif- 
icantly to Go(r) due to bending energies. However, the channel walls induce 
an apparent larger stiffness as shown in Fig. [2j Thus, analyzing experimental 
data with the result presented in Ref. [H] would lead to significant larger 
values for the persistence length of the polymer. Finite values of the poten- 
tial strength E, i.e., of k c effectively suppress additionally wavevectors with 
k < k c in the radial distribution function, so that only modes with k w k c 
contribute significantly to the partition sum. 

In order to compare the explicit expression for the radial distribution G(r) 
given by Eq. to experimental data we need to relate the average projection 
length < L z > to the contour length L. In Gaussian approximation on finds 
< L z >= L/(l + a' 2 ) with the variance 

<?max 

o 

of the first derivative of the filament position / with respect to the lateral 
position z. The upper border of the integral a' 2 g max = 27r/£ where £ denotes 
the length scale where thermal fluctuations disappear, which is either the 
monomer length of the considered protein or the spatial resolution of the 
experimental setup as in our case. Evaluating data from microscopy experi- 
ments one does not measure the contour length L of the filament but only the 
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Figure 3: The analytical result for the radial distribution function Go(r) 
given by Eq. <^ (solid line) describes accurately the measurements of the 
end-to-end distance of the filaments by fluorescence microscopy (symbols, 
see Ref. [6]). We find consistently the persistence length L p « 15 ± 3[i for 
all measurements whereas a theory neglecting confinement (dotted lines with 
E — and k = 19, k = 25, k = 35) does not agree with the experimental 
data. Moreover, comparing the data and the theory without confinement 
(E = 0) would imply a dependence of the persistence length on the channel 
width D. 

projection Lu of L on the focal plane. Analogous to L z one finds in Gaussian 
approximation for the averages the relation < L z >=< L\\ > /(l + cr /2 /2) 
which is used in Eq. ([3]) when comparing experimental data with the the- 
oretical result. Notice, that this correction remedies the violation of the 
inextensibility constraint \t(s)\ = 1 within Gaussian approximation. 

In Fig. [3] we compare the theoretical expression ([3]) to experimental dis- 
tributions measured by fluorescence microscopy from a single actin filament 
in a microchannel. The data and description of the experimental setup can 
be found in Ref. [6]. The widths of the channels range from D = 1.5/im over 
D = 5.8/xm to almost infinite widths where the filament lengths are in general 
larger than 10/im. We fit the end-to-end distribution function using the 
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potential strength E and the bending rigidity k as free parameters and find 
quantitative agreement between the theoretical prediction for confinement 
and the experimental data. In all cases we find consistently the persistence 
length Lp « 15 ± 3//m in good agreement with Lp = 18//m from Ref. [6] 
considering the difficulties in measuring persistence lengths. It is important 
to notice that the agreement between theory and experiment for different 
channel widths is found for identical persistence length Lp. Without taking 
the effect of confinement into account, i.e., a finite value for k c in Eq. ([3]) one 
would obtain significant different values for Lp. Fig. [3j(a) shows also that 
an unconfined filament (E = 0) cannot be used to describe the data. Tak- 
ing the height of the peak would systematically overestimate the persistence 
length and lead to unrealistic values. The quality of the fits with finite E 
indicates that an effective parabolic potential is appropriate, if the strength 
E depends on the channel geometry. Here the channel geometry is set by a 
rectangular cross-section with width D\ — D and height D2 ~ 1.4/zm. In a 
final step we derive an analytic expression to map the channel geometry on 
a parabolic potential with E(Di, D 2 ) which can be used in Eq. ([3]) even for 
non-parabolic confinement of the filament. 

For hard wall confinement in a channel with rectangular cross-section the 
potential reads U^ hw \f) = U x (f x )U y (f y ) with = for -f < / < f 

and Ui(fi) = 00 otherwise. To do the mapping we employ a technique from 
Ref. [IB] which proofed to work well for two dimensional membranes and 
thin films [TU] . It considers the root mean square fluctuations amplitude 
a xy = (fxy) given by o~~ 2 = 2 3 ^ 2 /3k 1 ^ 4: E 3 ^ 4: for parabolic confinement. 

The central idea is to calculate a 2 separately in Fourier space for short 
wavelength fluctuations and in real space for long wavelengths, respectively, 
and assuming that both approaches lead to the same value. In real space 
the polymer is to this end devided into independently fluctuating segments 
of size £||. This self-consistent ansatz leads to the implicit equation 



0= / <#(r-<> ™" . (7) 
J —00 

for the fluctuation amplitude <J x ,y Here, £y = A = v/^f is the correlation 
length parallel to the channel given by the deflection length A. In case of 
the hard wall potential U^ hw ^ we arrive at the equation exp (—1/(16//)) = 
v /7r/Ierf (l/yTB/t) for the ratio // = cr 2 /D 2 . Because // ~ 0.063777 is con- 
stant, one finds the scaling a oc D analogous to Odijk's scaling relation 
A oc Lp^D 2 / 3 derived for spherical channels [T7j. Moreover, one can deter- 

I 4 8 

mine the prefactor analytically by applying 4(3ELp/j,3 Da = 1 yielding the 
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expression 

A = f/5 = 2Llal = 20 LfD^ (8) 
V E 

for the deflection length. Notice, that for a one dimensional filament the 
ratio \x = a 2 /D 2 pa 0.063777 is much smaller compared to the value /i ~ 0.2 
for a two dimensional membrane |18j . 

Computing E for the parameters of the measured data set shown in Fig. [3] we 
find a good agreement with the experimental data, which makes our explicit 
result an important tool to analyse biophysical in vitro studies where more 
complex constraints are relevant due to the biological environments. For in- 
stance, spatial constraints can not only be induced by solid channel walls but 
also by forces between parallel filaments within the channel. Repulsive forces 
lead to additional extension whereas attractive interactions may cause the 
formation of bundles of semiflexible filaments. Our analytical self-consistent 
approach seems to exhibit a discontinuous unbinding transitions of bundles 
recently found in Monte-Carlo simulations [20]. Also in steady shear flow the 
conformational dynamics of individual polymer molecules can be visualised 
and the probability distribution for the polymer extension can be measured 
[2T1 [22] . The explicit expression (J3j) for the radial distribution function may 
be applied to distinguish between extensions due to shear flow or due to 
spatial confinement in channels. Recently, the fluctuations of semiflexible 
polymers were visualised in a nematic liquid crystal where confinement due 
to rod-like colloidal particles leads to an elongation of the filament [23] . The 
observed coil-rod transition makes our approach applicable so that the ex- 
plicit expression for the radial distribution function may be used to describe 
the elongation of the filaments beyond the tangent correlation function used 
in Ref. [23] • We expect that our main result, the explicit expression for 
G(R), can be used to analyze experiments of more complex situations such 
as solutions of filaments even in inhomogeneous channel geometries which is 
crucial for the understanding of the cytoskeleton in vivo. 

References 

[1] A.R. Bausch and K. Kroy, Nature Physics 2, 231 (2006). 

[2] F. Gittes, B. Mickney, J. Nettleton, and J. Howard, J. Cell Bio 120, 923 
(1993). 

[3] A. Ott et al., Phys. Rev. E 48, R1642 (1993). 

[4] J. Kas, H. Strey, and E. Sackmann, Nature 368, 226 (1994). 



9 



[5] L. Goff, O. Hallatschek, E. Frey, and F. Amblard, Phys. Rev. Lett. 89, 
258101 (2002). 

[6] S. Kdster, D. Steinhauser, and T. Pfohl, J. Phys. Cond. Mat. 17, S4091 
(2005). 

[7] T. W. Houseal, C. Bustamente, R. F. Stump, M. F. Maestre, Biophys. 
J. 56, 507 (1989). 

[8] M. Matsumoto et al, J. Polym. Sci. B Polum. Phys. 30, 779 (1992). 

[9] T. T. Perkins, D. E. Smith, S. Chu, Science 264, 819 (1994). 

[10] C. Bustamante, Z. Bryant and S. B. Smith, Nature 421, 423 (2003). 

[11] J. Han and H. G. Craighead, Science 288, 1027 (2000). 

[12] D. Nykypanchuk, H. H. Strey, and D. A. Hoagland, Science 297, 987 
(2002). 

[13] J. O. Tegenfeldt, et al. PNAS 101, 10979 (2004). 

[14] J. Wilhelm and E. Frey, Phys. Rev. Lett. 77, 2581 (1996). 

[15] M. Gardel et al., Phys. Rev. Lett. 91, 158302 (2003). 

[16] O. Kratky and G. Porod, Reel. Trav. Chim. Pays-Bas Belg. 68, 1106- 
1123 (1949). 

[17] T. Odijk, Macromolecules 16, 1340 (1983). 

[18] K. Mecke, T. Charitat, and F. Graner, Langmuir 19, 2080 (2003). 

[19] J. Vorberg, St. Herminghaus, and K. Mecke, Phys. Rev. Lett. 87, 196105 
(2001). 

[20] J. Kierfeld, T. Kiihne, and R. Lipowsky, Phys. Rev. Lett. 95, 038102 
(2005). 

[21] D. Smith, H. Babcock, S. Chu, Science 283, 1725 (1999). 

[22] C. M. Schroeder, H. P. Babcock, E. S. G. Shaqfeh, and S. Chu, Science 
301, 1515 (2003). 

[23] Z. Dogic et al., Phys. Rev. Lett. 92, 125503 (2004). 



10 



